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1.  INTRODUCTION 

This  report  contains  experimental  and  analytical  descriptions  of 
IEMP  produced  with  electron  beams.  The  experimental  work  was  done  for 
the  Harry  Diamond  Laboratories  under  contract  DAAG39-74-C-0006 . The 
data  reduction  and  analysis  were  performed  for  the  Defense  Nuclear 
Agency  under  contract  DNA001-75-C-0071 . 

The  basic  objective  of  the  work  was  to  obtain  experimental  data 
on  space-charge  neutralization  in  cavities,  to  test  the  adequacy  of 
the  existing  theoretical  model. 

The  Simulation  Physics,  Inc.,  SPI-5000  electron-beam  machine  was 
used  on  this  program,  modified  by  SPI  personnel  with  the  addition  of  a 
30 -cm- diameter  field  emission  diode  and  a vacuum-tight  transmission 
anode.  The  design  and  construction  of  the  test  chamber  and  the  inter- 
nal sensors,  including  a simple  magnetic  spectrometer,  were  performed 
at  IRT. 

Two  preliminary  experiments  were  conducted  using  4-  to  6-inch 
diameter  electron  beams  from  the  Simulation  Physics  SPI  Pulse  5000 
machine.  Results  are  summarized  in  References  1 and  2.  Rosado  (Ref. 
1)  measured  transmitted  beams  as  a function  of  chamber  depth.  Stahl 
(Ref.  2)  measured  transmitted  electrons  as  a function  of  chamber 
depth,  gas  pressure,  gas  composition,  and  collector  impedance.  An 
attempt  was  made  to  compare  the  measured  pulse  shapes  with  analytical 
results;  this  comparison  was  inconclusive. 

1.  J.  A.  Rosado,  "Space-Charge-Limited  Currents  in  Cavities," 

DNA  I CMP  Symposium  (U) , DNA  5098P,  June  1975  (SRD) . 

2.  R.  H.  Stahl  et  al.,  "Pressure  Effects  on  Space-Charge- 
Limited  Current  Transmitted  Across  a Cylindrical  Cavity,"  IE EE 
Trans.  Nucl.  Sci.  NS- 20 , No.  6,  December  1973. 
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The  experimental  effort  has  amassed  15  data  sets,  involving  some 
185  shots,  and  approximately  1000  oscilloscope  photographs.  From  this 
group,  four  selected  sets  (cases)  of  data  have  been  reduced.  In  each 
data  set,  the  cavity  geometry  and  parameters  of  the  injected  current 
are  held  constant;  the  varied  parameter  is  the  cavity  air  pressure, 


from  10"4  torr  or  lower  to  at  least  as  high  as  0.3  torr.  Table  1 sum- 


marizes 

the  four 

cases . 

Table  1 

SUMMARY  OF  DATA 

Case 

Number 

Shot 

Dates 

Shot 
Numb  ers 

Cavity 

Diameter 

Cavity 

Depth 

Emitted 

Current* 

FWHM 

Peak 

13 

8/74 

5178-96 

30 

cm 

14  cm 

30  nsec 

4 kA 

6 

8/74 

5235-42 

30 

cm 

30  cm 

30  nsec 

4 kA 

4 

8/74 

5214-25 

45 

cm 

30  cm 

30  nsec 

4 kA 

15 

I 

8/74 

5138-64 

30 

cm 

14  cm 

100  nsec 

10  kA 

Pleasured  with  1^, 

monitor. 

The 

cases  above  were  chosen  so 

that  three  comparisons  could 

be  made: 

1.  The 

effect  of 

cavity  depth  - case  13 

versus  case  6 

* 

■ 

2.  The 

effect  of 

cavity 

diameter 

- case 

6 versus  case 

4, 

3.  The 

effect  of 

emitted 

current 

- case 

13  versus  case  15. 

The  SPI 

field  emission  diode 

gap  was  1.5  cm 

and  the 

machine  charging 

voltage  was  250  kV  for  all  cases. 

The  experimental  results  are  now  in  a form  suitable  for  theoretical 
analysis,  and  a portion  of  the  analysis  is  included.  The  general  trends 
in  the  transmitted  current,  such  as  variation  of  the  pulse  shape  with  pres- 
sure, are  described  by  computer  simulation.  However,  some  of  the  details 
of  the  spatial  distribution  of  the  transmitted  current  are  not  presently 
understood . 


6 


2.  VACUUM  CHAMBER  AND  EXPERIMENTAL  CAVITY 


The  IEMP  cavity  is  an  adaptable  affair  that  is  contained  within  a 
larger  vacuum  chamber,  as  shown  in  Figure  1.  The  vacuum  chamber  attaches 
to  the  electron-beam  machine  via  an  insulating  ring  such  that  the  Mylar/ 
mesh  anode  of  the  machine  diode  becomes  the  emitting  end  of  the  cavity.  A 
removable  liner  and  a rear  collector  assembly  form  the  walls  and  rear  end 
of  the  cavity. 

The  liner  is  an  aluminum  cylinder  that  is  supported  and  electrically 
connected  to  the  chamber  walls  via  a collar  containing  rf  fingers  under 
compression.  The  rear  collector  assembly  consists  of  a central  circular 
aluminum  plate,  two  concentric  aluminum  rings,  and  a supporting  sleeve. 
Adjacent  collectors  are  supported  and  electrically  interconnected  by  a 
circular  ring  of  resistors.  The  support  sleeve  maintains  electrical  contact 
with  the  liner  through  rf  fingers.  The  depth  of  the  IEMP  cavity  can  be 
changed  simply  by  sliding  the  rear  collector  assembly.  The  liner  is  removed 
and  replaced  with  a short  liner,  which  is  moved  to  the  rear  to  create  the 
larger  diameter  cavity.  For  these  four  cases,  the  edge  of  the  liner  was 
positioned  1 cm  from  the  anode. 

The  vacuum  system  (not  shown)  consists  of  a 4-inch-diameter  oil  dif- 
fusion pump  and  roughing  pump.  The  cavity  pressure  was  monitored  by 
discharge  gages  and  a thermistor  gage,  which  have  been  calibrated  against 
a McLeod  pressure  gage. 

A sketch  of  the  magnetic  spectrometer  used  to  measure  electron  energy 
distributions  is  shown  in  Figure  2.  The  spectrometer  consisted  of  two  steel 


pole  pieces,  each  with  a 2.5-inch  radius  and  0.5-inch  thickness,  separated 
by  a 0.5-inch  gap.  A static  magnetic  field  in  the  gap  was  obtained  from 
several  permanent  magnets  attached  to  the  pole  pieces  across  the  gap  at 
the  outer  radius.  The  field  strength  in  the  gap  depended  on  the  number  of 
magnets  used.  The  ten  aluminum  collectors,  0.13  inch  wide  and  0.38  inch 
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to  a miniature  50-ohm  coaxial  cable  that  terminated  in  a BNC  feedthrough 
connector  at  the  wall  of  the  chamber,  where  the  oscilloscope  connection 
could  be  made.  The  spectrometer  assembly  was  attached  behind  the  rear 
collector  assembly,  with  its  0.12-inch  by  0.5-inch  rear  entrance  slit 
positioned  on  the  cavity  centerline.  All  measurements  were  made  at  high 
vacuum,  with  the  rear  collector  1 cm  away  from  the  Mylar  window. 

Table  2 gives  the  acceptance  range  of  electron  energy  at  each  col- 
lector for  the  two  static  field  strengths  used,  and  estimated  uncertain- 
ties . 

It  was  found  necessary  during  the  experiments  to  use  a relatively 
large  0.75-inch-diameter  front  aperture  (as  shown  in  Figure  2)  to  obtain 
reasonable  signal-to- noise  ratios.  This  degrades  the  energy  resolution 
(AE)  at  collector  #10  to  about  ±10%,  and  at  the  smaller  radii,  proportion- 
ately less.*  At  collector  #1,  the  effects  of  the  fringing  field  and 
space  charge  limiting  produce  an  estimated  ±10%  energy  resolution.  The 
results  of  the  measurements  are  discussed  in  Sections  3.1  and  3.4. 


Table  2 

SPECTROMETER  ELECTRON  ENERGIES 


Collector 

Number 

Collector 

Radius 

(cm) 

Electron 

Energy  1 

(keV) 

B 

gap 

= 

233  ±3  gauss 

B 

gap 

= 296  ±5 

gauss 

1 

1, 

.4 

7. 

.9- 

10. 

.0 

±10% 

12 

.7- 

15 

.9 

±10% 

2 

2, 

.0 

16 

.5- 

19 

.3 

17 

.8- 

30 

.8 

3 

2. 

.4 

25, 

.1- 

28 

.6 

40 

.0- 

45. 

.4 

4 

2, 

.8 

33 

.7- 

37 

.6 

53, 

.4- 

59 

.6 

5 

3, 

.1 

42. 

,7- 

47, 

.0 

67. 

.3- 

74, 

.0 

6 

3. 

.4 

51 . 

.3- 

56. 

.1 

80, 

.7- 

87, 

.9 

7 

3 , 

.7 

60, 

.0- 

65 

.2 

93 

.9- 

102 

.0 

8 

4 - 

.0 

68, 

.7- 

74  , 

.2 

107 

.0- 

115, 

.0 

9 

4, 

.3 

77, 

.8- 

83 

.6 

121 

.0- 

130 

.0 

10 

4. 

.5 

86, 

.8- 

92 

.7 

±10% 

134 

.0- 

143. 

.0 

±10% 

The  ±10"  energy  resolution  is  obtained  even  with  the  relatively 
wide  0.75- inch  aperture,  due  to  180°  focusing  inherent  in  the  design. 


3.  DATA  ORGANIZATION 


Table  3 is  an  index  of  the  eight  figures  relating  to  each  data  case. 
Note  that  Figures  C,  7,  and  8 describe  the  emitted  current  parameters  for 
all  three  short -pulse  cases.  The  measured  quantities  are  defined  in  the 
discussion  of  the  case  13  data,  Section  3.1. 


INDEX  OF  THE  FIGURES 


Figure  Titles 

Figure 

Parameters  Case  13 

Numbers 
6 4 

15 

Peak  currents  in  rear  and  wall 
collectors  versus  pressure 

IT,  IpC1>  IpC2> 
IFC3 

10 

14 

18 

Time  of  occurrence  of  peak 
currents  in  rear  and  wall 
collectors  versus  pressure 

II  ^ 

T’  FC1 

11 

15 

19 

Peak  current  density  averaged 
within  each  ring  of  rear 
collector  versus  pressure 

Jl’  J2’  J3  5~ 

12 

16 

20 

Current  collected  at  1 cm  depth 
(IpC1)  and  total  current  (IT) 

I^,  IpC1  at  1 cm  (i 

depth 

0 

6 

21 

Estimation  of  peak  current 
uniformity  through  trans- 
mission anode 

Jj,  J , J3  at  1 cm 
depth 

*7 

/ 

22 

Estimation  of  average  electron 
energies  versus  time 

V^,  IRT  Spectrometer  8 
data 

8 

8 

23 

Comparison  of  Ip^  and  1^ 
time  histories 

I I t) 

Fcr  T 

13 

17 

21 

3.1  CASE  13  DATA 

Figure  3 shows  the  peak  currents  collected  on  the  rear  and  the  wall  of 
the  test  cavity,  as  measured  in  the  four  shunt  resistor  rings  shown  in 
Figure  1.  The  currents  are  defined  as  follows. 


FC1 


the  total  current  collected  on  the  rear  plate. 


XFC2  and  TFC 3: 


the  currents  collected  within  a radius  of  10.2  cm 
and  5.1  cm,  respectively. 

I~:  the  current  collected  on  the  liner  plus  the  rear 

plate  current. 

The  shunt  resistances  are  each  approximately  10  mfi;  thus,  the  maximum 
voltage  drop  across  each  shunt  is  several  hundred  volts.  The  voltage  drop 
across  the  inductance  formed  by  the  liner  collar  and  the  1^  shunt  is  approxi- 
mately 600  volts.  The  shunts  were  calibrated  against  the  machine  current 
monitor  using  a resistive  load  between  the  cathode  and  the  center  of  the 
rear  collector,  and  passing  the  machine  current  through  all  shunts  in  series. 

In  August  1974,  a series  of  open-shutter  photographs  (reproduced  in  Appen 
dix  A)  was  made  of  the  beam  and  the  rear  collector  through  a chamber  window. 
Some  photographs  show  bright  spots  at  the  steel  screws  in  the  collector  rings 
and,  in  several  cases,  also  at  the  gap  in  front  of  the  Ip^  shunt.  The 
largest  spots  are  observed  at  the  100-nsec  pulse  width.  There  was  specula- 
tion that  the  Ic  _ measurement  was  disturbed  by  this  effect.  A review  of 
the  Ip£3  waveforms  recorded  during  this  series  do  not  show  any  unusual 
effects  when  the  local  optical  emission  is  observed.  It  is  possible  that 
the  optical  emission  is  associated  with  reocmbination  occurring  at  late  times. 

For  case  13,  the  pressures  were  varied  from  5 x 10  ^ to  20  torr,  in 
18  steps.  A majority  of  shots  was  confined  to  the  pressures  between  20  and 
150  mtorr.  Note  the  large  increase  in  currents  between  20  and  60  mtorr. 

See  Figures  9a-9d  for  the  time  histories  of  those  currents  indicated  by  a 
circled  data  point. 

Figure  4 shows  the  time  of  peak  currents  in  the  rear  collector  versus 
pressure.  Note  that  the  IT  and  Ip^  current  peaks  do  not  occur  at  the 
same  relative  time  within  the  pulse  except  near  the  40-  to  150-mtorr  pres- 
sure region.  The  times  to  peak  are  referenced  to  the  initial  rise  of  the 
diode  voltage  waveform  (Figure  8)  . 
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Figure  5 shows  the  peak  current  density  within  each  ring  of  the  rear 
collector  versus  pressure.  To  calculate  current  density  in  a ring,  the 
current  from  the  inner  rings  must  be  subtracted  before  dividing  by  the  ring 
area.  The  current  densities  J , J2>  and  J,  are  defined  below. 

t _ rFC3 

J5  2 

80  cm 


J2  = 


!FC2  ~ :FC3 
240  cm2 


TFC1  ' XFC2 
2 

400  cm 


Note  that  as  pressure  is  increased  beyond  45  mtorr,  the  current  densities 
J2  and  increase  significantly  more  than  . 

At  pressures  near  10  torr  and  above,  there  is  evidence  to  suggest  that 
the  Mylar/mesh  anode  may  have  bowed  slightly,  changing  the  emitted  current 
profile.  Both  the  diode  voltage  and  current  waveforms  and  the  transmitted 
current  waveforms  suggest  that  the  field  emission  diode  gap  decreases, 
particularly  on  the  axis,  due  to  the  forces  on  the  experiment  side  of  the 
window  (22  pounds  at  10  torr) . This  would  account  for  the  observed  decrease 
in  the  diode  impedance  and  increase  in  the  current  density  on  axis. 

The  next  three  figures  describe  the  emitted  current  for  the  short- 
pulse  cases.  Figure  6 is  an  approximation  to  the  emitted  current  as  all 
measurements  are  necessarily  collected  current  measurements.  The  time 
histories  of  1^  and  Ip^  are  compared  in  this  figure,  when  the  rear 
collector  was  moved  to  1 cm  from  the  anode,  flush  with  the  edge  of  the 
liner.  A cavity  depth  of  1 cm  is  sufficiently  small  to  have  eliminated 
longitudinal  space-charge  effects  at  the  peak  of  the  pulse.  Similar 
measurements  made  at  2 cm  and  5 cm  cavity  depths  produced  nearly  the 
same  1^,  and  Ip^  peak  values,  indicating  little  limiting  at  the  peak. 
Following  the  peak,  limiting  seems  to  occur  rapidly,  as  discussed  fur- 
ther below.  The  FWHM's  of  IpC1  and  1^,  are  25  and  32  nsec,  respectively. 
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Figure  7 shows  the  uniformity  of  currents  through  the  transmission 
anode.  The  current  density  is  given  here  as  a function  of  radius,  as 
measured  at  the  rear  collector  positioned  1 cm  from  the  anode.  Although 
the  radial  resolution  is  coarse,  there  are  higher  collections  near  the 
edge  of  the  beam,  while  the  collections  at  the  center  and  middle  rings 
are  lower. * 

This  data  was  recorded  after  a new  cathode  was  installed.  The  current 
densities  in  the  center  of  the  beam  observed  using  the  old  cathode  were 
significantly  lower,  and  have  been  recorded  both  on  the  rear  collector  and 
in  an  array  of  CaF  TLD's. 

Note  that,  although  this  data  was  taken  with  a 20-nsec  FWHM  emitted 
pulse,  the  profile  is  believed  to  be  similar  to  the  30-nsec  FWHM  cases. 
However,  the  magnitudes  in  Figure  7 must  be  multiplied  by  1.4  to  correspond 
to  the  peak  magnitude  shown  in  Figure  6. 

Figure  8 is  an  estimation  of  electron  energies  in  the  emitted  pulse. 

A plot  of  the  diode  voltage  waveform  is  given,  recalibrated  using  the  latest 
SPI  calibration  of  9.8  kV/V  and  corrected  for  a 22-nH  diode  inductance. 

This  curve  corresponds  to  electrons  impinging  on  the  Mylar  window.  A 
second  curve,  corresponding  to  electrons  leaving  the  window  and  entering 
the  test  cavity,  was  obtained  by  subtracting  the  electron  energy  losses 
in  the  0.25-mil  Mylar  window.  The  methods  used  to  calculate  energy  loss 
and  energy  spread  (AF.)  are  discussed  in  Section  4. 

Also  plotted  on  Figure  8 are  data  from  our  magnetic  spectrometer, 
showing  the  time  of  occurrence  of  the  peak  signal  versus  the  width  of  the 
spectrometer  energy  bln.  The  good  agreement  obtained  in  two  of  the  three 
bins  may  be  somewhat  fortuitous,  considering  the  poor  s i gnal -to-noi se 
ratios  of  the  short-pulse  spectrometer  oscilloscope  records.  However, 
since  much  better  agreement  between  spectrometer  and  modified  diode  wave- 
forms was  obtained  in  the  case  of  the  long  pulse  using  good-quality  records 
(see  Section  3.4),  we  tend  to  have  faith  in  the  modified  diode  waveform 
shown  here. 


Filter  measurements  by  SPI  personnel  confirmed  this  radial  distri- 
but ion  measurement. 


A comparison  of  Figures  6 and  8 shows  that  the  average  electron 
energy  at  50  nsec,  the  peak  of  the  I and  Ip^  pulses,  is  ^30  keV,  a 
value  high  enough  to  produce  little  limiting  at  1 cm  depth.  However, 
between  50  and  60  nsec,  the  average  energy  decreases  rapidly  to  several 
keV.  Ip^  decreases  rapidly,  while  the  wall  current  1^,  decreases  more 
slowly,  implying  a shift  from  longitudinal  to  transverse  electron 
trajectories. 

In  Figures  9a-9d,  the  time  histories  of  current  Ip^  and  1^  are 
plotted  for  seven  pressures.  The  oscillations  on  the  Ip^  waveform 
are  most  likely  a consequence  of  a small  resonant  network  between  the 
IpCP  shunt  resistor  ring  and  the  RF  finger  contacts  and  can  be  ignored. 
The  values  of  L and  C calculated  from  the  geometry  of  the  network 
yield  a resonant  frequency  of  280  mHz,  very  close  to  the  observed  ring- 
ing at  290  mHz.  Note  that  a small  prepulse  is  observed  at  7 to  8 nsec, 
very  likely  associated  with  the  initial  peak  of  the  diode  voltage.  The 
Mylar  anode  is  positively  charged,  and  is  the  source  of  displacement 
current  through  the  Faraday  cup  resistors.  The  main  current  pulse  rises 
'vlO  nsec  after  the  initial  rise  of  the  diode  voltage. 

3.2  THE  EFFECT  OF  CAVITY  DEPTH:  CASE  6 VERSUS  CASE  13 

The  case  6 initial  conditions  are  the  same  as  those  of  case  13, 
except  that  the  cavity  depth  is  increased  from  14  to  30  cm.  The  larger 
volume  thus  formed  might  be  expected  to  give  rise  to  a greater  negative 
potential  on  axis,  and  an  increased  electron  divergence  due  to  space- 
charge  buildup.  Comparing  Figure  10  to  Figure  3,  we  see  at  vacuum 
approximately  the  same  total  current  (1^)  but  significantly  less  rear 
plate  current  (Ip^),  as  expected.  As  pressure  increases  to  20  to  80 
mtorr,  self-neutralization  increases  for  both  cylinder  depths  and  a 
larger  percentage  of  the  charge  is  collected  at  the  rear.  It  is  inter- 
esting to  note  that  the  total  current  1^.  also  increases  in  this  pressure 
range.  It  seems  plausible  that  low-energy  electrons,  able  to  penetrate 
farther  into  the  cavity,  produce  additional  fields  which  cause  more 
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ionization  electrons  near  the  wall  to  accelerate  to  the  walls  and  be 
counted  in  1^.  At  the  higher  pressures  near  0.3  torr,  1^  decreases 
and  the  currents  1^  and  Ip^  are  very  similar  in  both  cases;  that  is, 
increasing  the  cavity  depth  has  not  significantly  reduced  the  total 
rear  collector  current  at  0.3  torr. 

The  current  density  plots  (Figures  12  and  5)  are  useful  in  examin- 
ing the  beam  current  profile.  Note  that  in  both  cases  the  collected 
beams  are  hollow  (J^  > ^ > jp  from  vacuum  to  ^50  mtorr.  As  pressure 
increases  above  50  mtorr,  both  plots  show  an  increase  in  ^ and 
(the  central  beam  density)  and  a lesser  change  in  (the  beam  edge 
density).  At  higher  pressure,  the  beam  in  case  6 (30  cm  depth)  appears 
to  be  focused  (J^  > J2  > jp  . The  profile  in  case  13  (14  cm  depth) 
to  5 to  6 torr  is  an  annulus  (Jj  > J3  > jp  , however.  It  has  been  sug- 
gested that  the  rate  of  pinch  is  such  that,  in  the  short  (14  cm)  cavity, 
the  beam  has  not  had  time  (or  distance)  to  pinch  to  less  than  the  FC2 
collector  radius,  while  in  the  30-cm-deep  cavity,  it  has  pinch  at  least 
to  the  radius  of  the  FC3  collector. 

A comparison  of  the  time  to  peak  current  data  of  Figures  11  and  4 
reveals  that  the  currents  IT  and  IpC1  peak  at  the  same  time,  only  near 
those  pressures  where  the  currents  are  maximum.  Otherwise,  IfC1  always 
peaks  before  IT  reaches  its  maximum  value*.  We  also  see  that  the  IpC1 
peak  occurs  much  earlier  in  the  30-cm  case.  This  is  a consequence  of 
the  onset  of  space-charge- limiting  earlier  in  time;  compare  the  vacuum 
time  histories  of  Figure  9a  and  13a. 

The  high-pressure  time  histories  of  Figure  13b  are  similar  to  those 
of  Figures  9b  and  9c,  except  that  1^  waveforms  show  significantly  longer 
decay  times  at  late  times.  It  may  be  that  these  "tails"  are  due  to  mag- 
netic effects  in  the  high-conductivity  plasma. 


*At  low  pressures  (pre-breakdown),  this  effect  is  due  to  space- 
charge-limiting.  At  high  pressures  (post-breakdown),  the  mechanism 
has  not  been  identified,  but  magnetic  field  (inductive)  effects  are 
suspected . 
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3.3  CASE  4 


Case  4 (Figures  14-17)  is  the  same  as  case  6 except  that  the 
cavity  diameter  has  been  increased  to  approximately  45  cm.  This  feat 
was  accomplished  by  removing  the  30-cm-diameter  liner  and  substituting 
a short  liner  on  the  collar  to  support  the  rear  collector  (see  Figure  1) . 
The  rim  of  the  short  liner  was  positioned  flush  with  the  surface  of  the 
rear  collector,  and  the  collar  slid  to  the  rear  of  the  chamber  to  in- 
crease the  cavity  depth.  Note  that  the  radii  of  the  emitter  and  the 
rear  collector  remain  the  same;  the  significant  change  is  the  removal 
of  the  conducting  cylinder  boundary  conditions  at  the  15-cm  radius. 

3.4  CASE  15 

Case  15  (Figures  18-24)  is  the  same  as  case  15  except  for  the  param- 
eters of  the  emitted  current  pulse.  Figures  21  and  22  characterize  the 
emitted  current  by  the  collected  current  in  a 1-cm-deep  cavity.  This 
pulse  is  significantly  wider  (100-nsec  BVHM)  and  more  intense  than  the 
short  pulse.  However,  the  I waveform  in  Figure  21  exhibits  the  same 
long  tail  observed  in  Figure  6 for  the  short  pulse,  implying  that  some 
low-energy  electrons  are  "squirted"  to  the  sides  of  this  1-cm-deep  cavity. 
The  current  density  versus  radius  plot  of  Figure  22  is  interesting,  show- 
ing a distinctly  different  beam  profile  from  the  short-pulse  cases.  Com- 
pared to  Figure  7,  the  long  pulse  exhibits  a higher  density  in  the  center 
and  an  improved  overall  uniformity. 

Figure  23  shows  electron  energies  as  a function  of  time  for  the 
long  pulse.  The  two  continuous  curves  are  obtained,  as  described  in 
Section  3.1,  from  the  diode  voltage  waveform  of  a typical  shot.  The 
data  points  are  obtained  from  the  magnetic  spectrometer  measurements 
representing  the  time  of  occurrence  of  the  peak  signal  placed  at  the 
bin  midpoint  energy.  The  agreement  is  very  good  except  at  very  early 
times.*  The  uncertainty  in  the  timing  accuracy  is  about  ±2  nsec. 

■k 

The  reason  for  the  modest  carly-time  discrepancy  is  uncertain; 
the  spectrometer  suffers  from  low  signal-to-noisc  ratios  at  early  times, 
whereas  the  voltage  probe  does  exhibit  some  overshoot  according  to 
recent  information  from  SP I. 
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The  uncertainty  in  the  midpoint  bin  energies  are  estimated  to  be  on  the 
order  of  ±10%,  as  discussed  in  Section  2.  Measured  values  of  AE 
( FWHM)  obtained  from  the  experiment  range  from  ±7.5  keV  at  45  keV  to 
±12.5  keV  at  100  keV. 


4.  DATA  ANALYSIS 


The  results  of  a computer  simulation  of  one  of  the  sets  of  data  are 
presented  here  and  compared  to  the  experimental  measurement.  Some  agree- 
ment is  obtained  for  the  current  at  the  outer  radius  of  the  collecting 
wall.  Time  histories  are  similar  and  trends  in  the  decay  of  the  persisting 
currents  are  duplicated.  However,  we  note  significant  differences  in  the 
vacuum  response,  as  well  as  in  the  spatial  distribution  of  peak  currents 
on  the  collecting  surface.  Whereas  the  differences  in  vacuum  response 
may  be  attributable  to  inaccuracies  in  the  electron  spectra  or  zoning 
used  in  the  calculation,  it  does  not  appear  that  the  difference  in  spa- 
tial distribution  can  be  easily  explained. 

4.1  APPROXIMATIONS  IN  THE  ANALYSIS 

The  data  of  case  13  were  analyzed  to  determine  to  what  extent  the 
effects  of  the  presence  of  the  background  gas  can  be  simulated.  The 
calculations  were  performed  using  the  DYNACYL  code  (Ref.  3).  In  light 
of  the  fact  that  the  method  used  to  treat  the  secondary  electrons  is 
strictly  applicable  only  when  the  concept  of  an  electron  drift  velocity 
is  meaningful,  and  the  fact  that  the  high  fields  encountered  here  cause 
runaway  so  that  the  concept  of  a drift  velocity  is  not  applicable,  a minor 
modification  had  to  be  made  to  the  code  to  allow  a meaningful  analysis  of 
the  experiments.  To  put  this  in  perspective,  a brief  summary  of  the  treat- 
ment is  given  here. 

In  these  calculations,  the  secondary  electrons  are  treated  as  a fluid. 
Primary  electrons  cause  ionization  of  the  gas  at  a rate  determined  by  the 


3.  E.  P.  dePlotnb  et  al.,  "Two-Dimensional  Time-Dependent  Computer 
Code  Development  and  Applications  to  Problems  in  IEMP  (U),"  DNA  35051-', 
August  1974  (SRD) . 
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Bethe  approximation  for  slowing  down,  and  the  conversion  of  primary  energy 
to  ionization  is  at  the  rate  of  three  Rydbergs  per  ion  pair.  This  has  the 
limitation  that  the  secondaries  dissipate  their  initial  high  energy  in 
producing  further  ionization,  possibly  not  the  case  for  extremely  high 
ratios  of  electric  field  to  gas  number  density  (E/N).  In  the  latter 
case,  this  ionization  rate  may  be  high  by  as  much  as  a factor  of  two. 

The  fluid  of  secondaries  is  followed  in  an  Eulerian  grid  with  drift 
velocity  determined  by  experimental  values  of  the  velocity  of  electrons 
in  nitrogen.  This  assumes  that  the  collision  frequency  for  momentum 
exchange  is  large  compared  with  other  frequencies  of  interest.  As  to  fur- 
ther ionization  by  secondaries  (the  cause  of  avalanche) , the  treatment  is 
again  based  on  experimental  data.  What  is  more,  this  empirical  treatment 
shows  a phenomenon  which  is  the  basis  for  the  method  used  to  extend  the 
calculations  to  the  regime  covered  by  the  present  experiments.  This  phe- 
nomenon is  the  decrease  of  ionization  rate  above  a certain  value  of  E/N. 

The  secondary  ionization  rate  shows  an  increase  with  E/N  up  to  10*7 

2 

V-m  . At  that  point,  the  rate  begins  to  decrease  with  increasing  field. 
This  is  at  a point  where  the  drift  velocity  begins  to  exceed  10  m/sec. 

To  maintain  the  simplicity  of  the  method,  an  ad  hoc  value  for  the  drift 
velocity  was  chosen  for  fields  causing  higher  secondary  velocities.  In 
particular,  we  choose  a drift  velocity 

(v  + 0.5c  x\  . 

— — : 1 for  v > 10  m/sec  , 

1 + x / e 

where  v&  is  the  tabulated  velocity,  c is  the  velocity  of  light,  and  the 
dimensionless  quantity  x is 

x = 1017  E/N  . 

The  direction  of  drift  is  parallel  to  the  electric  field.  Although  this 
ignores  the  inertial  properties  of  the  secondaries,  it  maintains  an  essen- 
tial physical  characteristic  of  the  problem:  the  secondaries  travel  rapidly 
in  the  direction  of  the  fields.  However,  since  the  effects  of  magnetic 
forces  are  ignored,  the  assumption  of  motion  parallel  to  the  electric 
field  introduces  a source  of  error. 
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There  are  other  assumptions  and  simplifications  incorporated  into 
the  calculation  (Ref.  3) , but  the  use  of  a drift  velocity  is  the  most 
suspect  in  the  present  analysis. 

4.2  RESULTS  OF  THE  ANALYSIS 

Seven  experiments  were  analyzed.  The  experimental  results  are  pre- 
sented in  Figures  9a-9d.  Calculations  were  performed  for  a tank  15  cm  in 
length  and  15  cm  in  radius.  The  problem  was  represented  by  a grid  of  five 
zones  in  length  and  five  in  radius.  This  limited  zoning  was  imposed  by 
the  restriction  of  maintaining  a tractable  quantity  of  output  and  the 
reluctance  of  the  analyst  to  permit  particles  representing  electrons  to 
traverse  more  than  a zone  in  a particle  time  step.  However,  the  coarse 
zoning  may  cause  substantial  error  in  the  vacuum  and  20-mtorr  cases, 
and  at  higher  pressures  for  early  times  (prior  to  the  peak  of  the  pulse). 
Emission  took  place  from  45  equally  spaced  radial  points  on  the  emission 
surface  for  the  two  lowest  pressures  and  15  points  for  the  five  higher- 
pressure  cases. 

For  emission  from  zones  within  10.5  cm  of  the  axis,  the  current  was 

reduced  by  a factor  of  0.3  relative  to  the  emission  at  radii  greater  than 

10.5  cm.  This  was  motivated  by  the  curve  of  measured  currents  shown  in 

2 

Figure  7.  The  peak  emission  current  density  was  7 amp/cm  . Current  time 
history  was  as  in  Figure  6 (Ip^) , and  the  digital  representation  is  shown 
in  Figure  25.  The  energy  time  history  was  taken  as  the  corrected  diode 
voltage  Vp  of  Figure  8,  modified  by 

v2‘ 


E = V. 


. /20  keV\2 

. ' \ VD  ) 


eV 


where  20  kcV  was  assumed  to  be  the  cutoff  of  the  Mylar  window.  This  for- 
mulation is  based  on  the  fact  that  electron  energy  loss  is  approximately 
inversely  proportional  to  the  electron  energy.  The  digitized  representa- 
tion is  shown  in  Figure  26. 

The  electron  emission  was  taken  to  be  uniform  in  solid  angle  confined 
to  the  polar  angle  interval  0 to  0 


max 


where  0 is  given  by 
max  6 1 


max 


_ H * 20  keV* 

= min  1.0,  — g , 

- — 
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where  0 is  in  radians.  This  is  an  approximation  to  a function  indi- 
cated by  SPI  and  is  shown  in  Figure  27;  we  note  that  the  spread  in  polar 
angle  appears  to  be  substantially  broader  than  what  one  would  expect  to 
see  from  electrons  incident  normal  to  the  Mylar  being  scattered  off  normal 
by  the  Mylar.  Furthermore,  some  studies  were  made  using  random  energy 
distributions;  these  will  be  discussed  after  the  main  body  of  results. 

The  results  for  the  seven  pressures  are  seen  in  Figures  28a-28c. 

The  calculated  results  are  at  13.5  cm  and  were  extended  to  15  cm  by  using 
linear  extrapolation  in  radius  by  the  formula 

I (r  = 15  cm)  = 2tt  15  H(r  = 13.5  cm) 

or 

I(r  ± 15  cm)  - H(r  = 13.5  cm)  x 1 meter  , 

where  I is  the  current  collected  on  the  rear  plate. 

Immediately  evident  from  Figure  28a  is  the  fact  that  predicted  results 
for  the  vacuum  case  are  higher  than  experiment  by  up  to  a factor  of  two. 
This  could  be  explained  if  the  true  electron  energy  were  lower  than  that 
used  in  the  calculations.  There  is  also  the  possibility  that  the  zoning 
is  too  crude  to  describe  the  two  lowest-pressure  cases.  In  the  vacuum 
and  20-mtorr  cases,  the  current  oscillated  over  the  range  indicated  in  the 
figure  and  has  been  averaged  over  4 nsec  for  presentation. 

At  a pressure  of  20  mtorr,  the  calculated  response  is  still  too  large, 
which  could  be  due  to  spectral  differences,  to  the  coarse  zoning,  or  to  the 
treatment  of  the  primary  ionization  rate.  The  last  possibility  is  examined 
later.  The  late-time  tail  is  represented  well  by  the  calculation,  a fact 
which  is  probably  fortuitous.  However,  the  shape  of  this  tail  is  probably 
valid  since,  after  collapse  of  the  fields,  the  drift  approximations  should 
not  be  so  bad. 

At  45  mtorr,  the  peaks  agree  well,  a fact  which  is  probably  not  too 
significant  since  nearly  all  of  the  current  is  getting  through  in  both 
experiment  and  theory.  The  tail  falls  off  a bit  too  rapidly  in  the  cal- 
culation, but  still  the  shapes  of  the  curves  are  similar. 
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The  90-mtorr  case  is  interesting  because  the  decrease  at  60  nsec  in 
the  experiment,  due  to  the  drop  in  driving  voltage,  is  not  predicted  by 
the  theory.  This  difference  can  be  explained  by  the  approximate  primary 
ionization  rate  being  high  (See  page  14).  Except  for  this  difference, 
the  agreement  is  good. 

Looking  at  the  150-mtorr  results,  it  can  be  seen  that  the  experimental 
results  also  lose  the  decrease  at  60  nsec,  a further  indication  that  the 
phenomena  are  being  described  accurately  but  begin  occurring  at  somewhat 
too  low  a pressure  (i.e.,  90  mtorr)  in  the  analysis.  The  theoretical 
falloff  in  the  tail  has  the  correct  slope. 

At  higher  pressure  (3.5  torr) , the  rate  of  decay  of  the  tail  has 
increased,  a fact  which  appears  in  both  theory  and  experiment. 

At  the  highest  pressure  examined,  15  torr,  the  agreement  deteriorates 
although  both  theory  and  experiment  show  a slower  decay  rate.  This  diver- 
gence is  not  considered  too  serious  since  the  high  pressure  on  the  anode 
may  deform  the  diode  and  change  the  spatial  distribution  of  the  electrons 
incident  on  the  chamber.  Since  the  calculated  tail  disappears  for  totally 
hollow  emission,  one  is  led  to  believe  that  the  falloff  would  be  slower 
in  the  case  of  more  uniform  emission,  a fact  which  is  consistent  with  the 
assumed  distortion  in  the  anode." 

In  the  calcuations  discussed  so  far,  we  considered  monoenergetic 
electrons,  with  electron  energy  varying  during  the  pulse.  To  calculate 
the  effect  of  an  energy  distribution  in  place  of  the  monoenergetic  elec- 
trons, we  performed  calculations  assuming  a distribution  approximated  by 


This  is  shown  in  Figure  29  along  with  the  data  supplied  by  SPI.  The 
energy  was  taken  as  random  about  the  central  energy  given  in  Figure  26 
with  a full  width  of  AE. 

For  the  vacuum  case,  the  oscillations  and  response  were  both 
decreased  by  approximately  20"o  by  the  use  of  an  energy  distribution; 
see  Figure  28a.  The  agreement  is  not  substantially  better.  For  the 
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20-mtorr  case,  the  ionization  rate  by  primaries  was  reduced  by  a factor 
of  two  to  examine  the  possibility  that  the  secondaries  do  not  imme- 
diately slow  by  further  ionization.  This  modification  and  the  use  of 
the  random  energy  distribution  produced  the  results  shown  in  Figure 
30,  which  may  be  construed  as  an  improvement  but  could  be  explained  by 
the  use  of  an  electron  spectrum  substantially  harder  than  the  real 
spectrum.  Again,  effects  of  the  coarse  zoning  may  be  significant. 

The  comparison  between  the  calculated  and  measured  current  on  the 
outside  of  the  collector  (Figures  28a-c)  is  not  too  discouraging.  How- 
ever, the  situation  deteriorates  when  one  examines  the  spatial  distri- 
bution of  current.  In  general,  the  current  calculated  theoretically 
shows  a definite  pinching  effect.  The  calculated  peak  current  enclosed 
with  a radius  of  4.5  cm  is  more  than  that  experimentally  measured  for 
the  5.2-cm  collector  (see  Figure  3).  As  an  example,  Figure  31  shows 
the  calculated  spatial  distribution  of  currents  for  the  150-mtorr  case, 
where  the  predicted  total-current  time  history  agreed  well  with  experi- 
ment. The  theoretical  values  lie  above  the  experimental  curve  for  the 
peak  interior  currents.  The  peak  currents  occur  at  successively  later 
times  for  decreasing  radii  which  may  be  due  to  pinching.  Since  magnetic 
forces  are  not  allowed  to  act  on  secondaries,*  there  was  interest  in 
studying  the  effect  on  the  calculations  of  ignoring  these  forces  on  the 
primary  particles  also.  To  investigate  this  possibility,  the  same  pro- 
blem (Figure  31)  was  run  ignoring  the  magnetic  forces  altogether.  As 
seen  from  Figure  32,  the  currents  maintain  their  hollow  shape.  Thus, 
it  appears  that  ignoring  the  effect  of  magnetic  force  on  the  secon- 
daries is  somewhat  compensated  for  by  the  neglect  of  this  force  on  the 
primaries;  but  there  is  no  longer  agreement  with  the  rate  of  decay  of 
the  currents. 

The  above  comparison  of  theory  and  experiment  leads  to  the  follow- 
ing conclusions.  The  method  used  here  of  treating  secondaries  as  moving 
at  a drift  velocity  outside  the  regime  of  strict  applicability  of  the 
concept  produces  satisfactory  results  for  total  transmitted  currents  at 

This  would  have  made  the  conductivity  a nondiagonal  tensor,  which 
the  code  could  not  handle. 
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pressures  above  45  mtorr  for  the  experiments  examined  here.  Detailed 
results  such  as  spatial  distribution  of  the  transmitted  current  are 
substantially  in  error.  In  particular,  a substantial  pinching  effect 
is  noticed  in  the  calculation  which  is  not  present  in  the  experimental 
measurements.  This  treatment,  which  was  initially  intended  to  comple- 
ment the  treatment  of  secondaries  as  inertial  particles  and  thereby 
eliminate  stability  problems  caused  by  plasma  oscillations,  does  not 
appear  adequate  to  describe  the  full  effects  of  the  background  gas. 
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Figure  I.  Vacuum  chamber  and  experimental  cavity 
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•igure  2.  IRT  magnetic  spectrometer 
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eurc  4.  Time  of  occurrence  of  peak  currents  in  rear  an 
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Current  collected  at  1 cm  depth  (Ipci)  and  total  current  ( 1 ) 
approximation  to  emitted  current  - short  pulse 
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f igure  8c.  Comparison  of  1^^.^  and  1^  time  histories 
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Peak  currents  in  rear  and  wall  collectors  versus  pressure 


Time  of  occurrence  of  peak  currents  in  rear  and  wall  collectors  versus  pressure 


current  density  averaged  within  each  ring  of  rear  collector  versu 
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Figure  17b.  Comparison  of  I„r,,  and  I...  time  histories 
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Figure  21.  Current  collected  at  1 cm. 
depth  ( T pQi ) and  total  current  (Iy); 
approximation  to  emitted  current-long  pulse 
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Tig ure  25.  Digitized  pulse  shape  used  to  describe  emission 

current  as  a function  of  time.  Th>  ' is  an  approx 
imation  to  the  curve  Ipr j of  Figure  6. 
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Figure  26. 


Digitized  energy  of  emitted  electrons  as  a function 
of  time.  This  is  an  estimate  of  the  electron  energy 
(Figure  8)  as  modified  after  passing  through  the 
Mylar  anode. 
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H (kA/m)  AT  OUTSIDE  OF  COLLECTOR 


Figure  30.  Effect  of  varying  the  ionization  cross  section. 

Reducing  the  ionization  reduces  the  transmitted 
current,  hut  the  injected  electron  energy  appears 
to  be  too  high. 
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Figure  31.  Spatial  distribution  of  collected  current  at  150 

mtorr.  This  demonstrates  that  the  measured  values 
are  much  less  in  the  center  than  the  calculated 
(solid  curves)  values.  Measured  values  are  from 
Figure  3 at  amplitudes  designated  by  the  arrows. 
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INCLUDED  CURRENT  ( kA ) 


Figure  32.  Spatial  distribution  of  collected  current  at  150  mtorr. 

Curves  are  calculated  values.  The  magnetic  forces  were 
ignored.  This,  in  comparison  with  Figure  3,  demonstrates 
the  pinching  of  the  magnetic  field.  The  detailed  shapes 
of  the  curves  for  radii  <13.5  cm  are  not  meaningful  since 
data  were  available  at  every  5 nsec  and  the  statistics 
are  worse  than  in  the  previous  figure  due  to  the  reduc- 
tion in  current.  Measured  values  are  from  Figure  3 at 
amplitudes  designated  by  the  arrows. 


APPENDIX  A 

OPEN -SHUTTER  PHOTOGRAPHS 


A series  of  open-shutter  photographs  was  taken  of  the  cavity  during 
the  e lectron-beara  measurements.  The  cavity  was  viewed  through  a window 
in  the  side  of  the  vacuum  tank.  A Polaroid  camera  fitted  with  a closeup 
lens  was  used.  The  illumination  is  due  solely  to  electron-beam  interac- 
tions. Thus,  the  photographs  give  a quantitative  time- integrated  picture 
of  the  progress  of  the  electrons  via  light  emitted  as  a result  of  colli- 
sions with  the  gas  and  cavity  materials. 

The  photos,  which  are  reproduced  on  the  following  pages,  were  taken 
during  a pressure  run  using  a long-pulse  input  and  a 30-cm-deep  cavity. 
The  cavity  diameter  was  45  cm  as  a result  of  removing  the  30-cm-diameter 
liner.  Pressures  range  from  14  mtorr  (Figure  A-la)  to  6 torr  (Figure 
A-lf).  At  14  mtorr,  the  rear  collector  and  side  wall  that  are  clearly 
illuminated  at  higher  pressures  are  obscured,  and  only  the  epoxy  surface 
of  the  magnetic  field  sensor  is  seen.  Here,  the  beam  has  "blown  up"  due 
to  space-charge-limiting.  At  32  and  60  mtorr,  the  beam  is  confined  to 
the  cavity  axis,  except  for  the  streamers  observed  in  mid-cavity.  At 
100  mtorr  and  above,  the  beam  appears  more  uniform  and  the  bright  emis- 
sion discussed  in  Section  3.1  is  observed  on  the  rear  collector  near  the 
screws  and  the  edges  of  the  collector  rings. 
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